WU-AP/51/96, |astro-ph/ 9601056 



Fluctuations of the gravitational constant 
induced by primordial bubbles 

Nobuyuki Sakai* 

Department of Physics, Waseda University, Shinjuku-ku, Tokyo 169, Japan 

(Revised April 10 1996) 

Abstract 

We consider the classical fluctuations of the gravitational constant generated 
by bubbles in the inflationary universe. For extended inflation, we demon- 
strate numerically how and how large fluctuations are produced during bubble 
expansion. The amplitude of the fluctuations depends on the Brans-Dicke pa- 
rameter u: if u> is of the order of unity, the amplitude becomes of the order 
of unity within one Hubble expansion time; if u is large (say, u = 1000), the 
growth rate of the fluctuations is small, but it keeps growing without freez- 
ing during inflation. We also discuss some astrophysical implications of our 
results. 
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Particle physics predicts that the universe experienced many phase transitions in its 
early history. If any of these phase transitions is first order, the universe changes its phase 
from false vacuum to true vacuum through the creation, expansion and collision of bubbles. 
Old inflation is based on a super-cooled first-order phase transition. In the inflation- 
ary scenario, the universe expands exponentially with time before the transition, thereby 
solving the horizon, flatness, and monopole problems. However, it turns out that this expo- 
nential expansion is too rapid to permit a transition from false vacuum to true vacuum via 
percolation of true vacuum bubbles ||]. 

Extended inflation || revived the idea of old inflation by using the Brans-Dicke the- 
ory instead of the Einstein theory. The Brans-Dicke field decelerates the expansion of the 
universe so that true vacuum bubbles can coalesce, thus ending the phase transition that 
drives inflation. This model provides an interesting hypothesis that the large-scale structure 
of galaxy distribution is generated by primordial bubbles [|J. Many discussions about ex- 
tended inflation have been made in past years. La et al. and Weinberg H estimated volume 
fraction of bubbles and found a difficulty: the constraint from the isotropy of the cosmic 
microwave background requires the Brans-Dicke parameter u < 25, which contradicts the 
lower limit u > 500 from the time-delay experiments ||. This problem is called a "big- 
bubble problem" because many big bubbles cause large-scale inhomogeneity. Thus interest 
shifted to other extended models based on general scalar-tensor theories [[/J- ||. 

In discussing density fluctuations in extended inflation or other extended models, several 
authors have estimated volume fraction of bubbles as we mentioned above, or quantum 
fluctuations of the gravitational constant (the Brans-Dicke field) jnj. In this paper, we 
investigate the classical fluctuations of the gravitational constant generated during bubble 
expansion. It was shown that, under the thin- wall approximation, the Brans-Dicke field 
must be inhomogeneous inside a bubble from the consistency of the junction conditions 



Hfl . However, it is not clear how and how large fluctuations are really generated. If these 



fluctuations are not small, we should consider it for the constraint for the models, such as 
the structure formation process and the measurements of the gravitational constant. 

There is an alternative scenario of inflation which is accompanied by bubble nucleation: 



one-bubble inflation |12| . This model is distinguishable from extended inflation or its gener- 
alized version: the second slow-rollover inflation occurs inside a nucleated bubble, and our 
observable universe is entirely contained in one bubble. If we assume the 0(4)-symmetric 



bubble, the interior of the bubble can be a homogeneous and isotropic open universe |13] ; this 
model may account for the universe model with f2 <0.1, which is supported by increasing 
observations. Although the fluctuations of the gravitational constant may not be relevant 
to this model, it is interesting to study the fluctuations inside a bubble in a similar way. 
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To begin with, we explain briefly why fluctuations of the gravitational constant are 
generated in the inflationary universe even at a classical level. Let us consider the system 
which is described by the action 



S = J dV^fjJ^STC - ^(V$) 2 } - U(9) - \{Vipf - Vty) 



(1) 



where $ is a Brans-Dicke-like scalar field which is normalized to be unity at the present 
epoch, and if) is a Higgs-like inflaton field. The field equation for $ is derived from ([!]): 
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2cu($) + 3 



-87rG {(WO 2 + V(if))} - u'{<S>){V<S>) 2 + £/'($)] . (2) 



Just looking at (0), we can understand that the inhomogeneity of if> causes that of <3>: 
even if $ and if) are homogeneous at the beginning of inflation, bubble nucleation, which is 
realized by the inhomogeneization of if), makes $ inhomogeneous inside a bubble. Because 
most bubbles become superhorizon scale at the end of inflation, the fluctuations of the 
gravitational constant, which trace these bubbles, cannot disappear soon after inflation. 

We may note, in passing, that the inhomogeneization inside a bubble occurs in generic 
systems. After a conformal transformation and the introduction of a new scalar field (f>, 



many scalar-tensor theories have the Lagrangian form |14 
S — I d 4 xy/—g 



(3) 



where j3 and 7 are dimensionless coupling constants, and k 2 = 8ttG. The field equation for 
is written as 

□0 = _Z|V 7 ^(V^) 2 - pKe-^V^) + W'(<p). (4) 

This equation also indicates that the inhomogeneity of if) causes that of 0. Thus we see that 
the inhomogeneity inside a bubble is a generic property in generalized Einstein theories or 
in multiple scalar field theories. 

In what follows we numerically investigate how and how large fluctuations of $ are really 
generated in the inflationary era. In order to solve the field equations, we must specify a 
gravitational theory, i.e., the functions of u;($) and U(<&). Here we simply assume u;($) = 
constant and U($>) = 0, which corresponds to extended inflation. Although our assumption 
is just for simplicity, our results will be also applied to some of the models with variable 
u, so-called hyperextended inflation []|||. Steinhardt and Accetta and Garcia-Bellido 
and Quiros Q independently proposed similar scenarios: a;($) is small and almost constant 
during inflation and diverges as $ approaches $o(= 1)- As long as we discuss their models, 
the assumption of constant u during inflation is justified. 
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We assume a spherically symmetric spacetime: 



ds 2 = -dt 2 + A 2 {t, r)dr 2 + B 2 {t, r)r 2 {d6 2 + sin 2 6d<p 2 ), (5) 



As a potential of the inflaton field, we adopt the quartic potential which was used in |15|: 



V(i>) = {ip- m) 2 {3(e + l)ip 2 + 2emip + em 2 } . (6) 

Free parameters are e and m: e determines the shape of the potential and the ratio of m to 
the Planck mass mpi implies the strength of gravity. We depict some shapes of the potential 
in Fig. 1. 

As an initial configuration of the inflaton field, tp(t = 0,r), we give a static bubble 
solution in a flat spacetime by numerically solving the field equation. As for the Brans- 
Dicke field, we assume $(i = 0, r) to be homogeneous. We also suppose that dip/dt(t = 
0, r) = d§/dt(t = 0,r) = 0. To solve the time-dependent field equations, we use a finite 
difference method with 1000 meshes; the details of the method were shown in the Appendix 



of ||16|| . The Hamiltonian constraint equation is used for checking the numerical accuracy; 
through all the calculations executed here, the errors are always less than a few percent. 

Our numerical results are summarized in Figs. 2—5. We normalize the time scale and 
the spatial scale by x~ l = {87rGV A (0)/3$(0)} _ 5, which corresponds to the horizon scale at 
the nucleation time. In Fig. 2 we show examples of the time-evolution of ip and $. We see 
that the fluctuations of <3> are really generated as a bubble expands. 

Here we define a quantity which corresponds to the amplitude of the fluctuations as 

_ $(t,r max ) -$(t,0) 
l j ~ $(t r ) ' { > 

Y V 1 1 max ) 

where r max is the outer numerical boundary. We utilize A(t) to present the following results. 
In Fig. 3(a) we draw A(t) for various lu. As we can expect from (0), A(t) at fixed t is almost 
proportional to 1/u: if u — 1, the amplitude becomes of the order of unity within one 
Hubble expansion time; if to = 1000, the growth rate of the fluctuations is small, but A(t) 
goes on increasing without freezing until one bubble collides with another. 



Samuel and Hiscock (15| investigated Euclidean solutions of an 0(4)-symmetric bubble 
and found their configurations for m/mpi of the order of unity are quite different from those 
of ordinary bubble solutions. Here we pay attention to the behavior of §(t, r) for large 
m/mp\. Figure 4 shows the evolution of $(t,r) for m/mp\ = 0.5. The behavior of $(t,r) 
is quite different from that for small m [cf. Fig. 2 (b)]. In Fig. 4(b) we draw A(t) for 
m/mpi = 0.1, 0.2, and 0.5. As m/mpi becomes larger, the fluctuations of $ become larger. 
We can interpret our results as follows. In the case of large m/mpi, the viscosity term of 
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d^/dt in the right-hand side of (|2|) is dominant to spatial derivatives, and therefore the 
homogeneization process is slowed. 

We also investigate the dependence of the fluctuations of $ on the potential shape. We 
draw A(t) for several e in Fig. 5. The potential for e = 0.05 corresponds to the case of a 
thin-wall bubble, and e = 0.5 to a thick-wall bubble. The temporal behavior of A(t) for 
small e (thin-wall bubble) is violent, but the eventual behavior depends little on e. 

So far we have numerically analyzed the evolution of a bubble in extended inflation 
and see how and how large fluctuations of the Brans-Dicke field are generated. In order 
to understand the astrophysical implications of our results, we also have to analyze their 
subsequent evolution after the phase transition. Because most fluctuations are nonlinear and 
gravity is not described by the Einstein theory, however, such analysis is not easy: another 
simulation or some sophisticated approaches are needed. In the rest of the paper, we offer 
discussion on the rough implications to some inflationary models. 

In the scenario of extended inflation or hyperextended inflation, all bubbles are nucleated 
with the subhorizon size and cross outside the horizon during inflation. After the phase 
transition, matter trace the bubbles— we call such bubblelike distribution after inflation 
"voids"— and reenter again during the radiation-dominated era or the matter-dominated 
era. Bubbles nucleated the earliest cross outside the horizon first, reenter last, and become 
the largest voids at present. Although we do not know the exact evolution of voids, we 
may suppose that homogeneization inside the void does not occur effectively until the void 
reenters the horizon. Because astrophysical-sized voids (30 — lOOMpc) reenters the horizon at 
z ~ 10 4 to 10 3 , the inhomogeneity of the gravitational constant affects the evolution of voids 
during the radiation-dominated universe. In relation to the big-bubble problem of extended 
inflation, Vadas [I?]] studied thermalization of superhorizon voids in the radiation-dominated 
universe and found homogeneization can occur in a much shorter time than previously 
thought. This result indicates that there remains a possibility that the big-bubble problem 
is resolved; it is worth studying further. In order to analyze the thermalization process 
more precisely, it is important to take the inhomogeneity of the gravitational constant into 
account. 

Let us discuss the following question: is it possible to observe the fluctuations of the 
gravitational constant at present? According to the linear perturbation analysis in the 
matter-dominated Brans-Dicke universe |18[ , the decay rate of the amplitude of $ is ap- 
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proximately proportional to (scale factor) - * for u ^> 1. This indicates that, if Scj)/(f> < 1 
when a void enters the horizon and the linear perturbation theory can be applied, the present 
fluctuations of <3? inside astrophysical-sized voids are too small to be observable. However, if 
the fluctuations are nonlinear at the time, it may be possible that a large inhomogeneity still 
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remains although the amplitude cannot be estimated without further analysis. Small-scale 
fluctuations caused by small bubbles are also interesting to study. While Turner and Wilczek 
p~9|] pointed out the detectability of the gravitational waves produced by bubble collisions 
in extended inflation, small-scale bubbles could be a source of "scalar-type" gravitational 
waves which can be detected by a laser interferometer. 

Finally, we make some comments on one-bubble inflation. If we consider a single scalar 
field in the Einstein theory, which is described by the 0(4)-symmetric bounce solution [ITS 



the interior of the bubble can be a homogeneous, isotropic, and open universe. We may 
understand that the high symmetry allows the inside spacetime to have a homogeneous 
and open slicing by a miracle. On the other hand, if we suppose two scalar fields or a 
modified Einstein gravity, 0(4) symmetry is lost and the inside must be inhomogeneous. 
Our analysis gives a generic picture that in such a theory the fluctuations grow during 
bubble expansion. Unless the amplitude of the fluctuations are negligibly small at the end 
of inflation, the inhomogeneity can be observed at present because it is still on a super- 
horizon scale. Recently, Linde and Mezhlumian [2(J proposed some models with two scalar 
fields, taking into account the condition that the inhomogeneity inside a bubble should be 
negligible. It may also be important to consider the classical fluctuations in models with a 



modified Einstein gravity such as proposed in |H] . 
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Figure Captions 

FIG. 1. The potential of the inflaton field we used in our analysis. We depict the shapes 
for e = 0.05,0.1, and 0.5. 

FIG. 2. Examples of the evolution of (a) ip and (b) $. We see how the fluctuations of $ 
are generated as a bubble expands. 

FIG. 3. Dependence of the fluctuations of $ on u. We draw the amplitude A(t) for 
w = 1, 10, 100, and 1000. We set e = 0.1 and m = 0.1. A(t) at fixed t is almost 
proportional to but it goes on increasing even for large uj. 

FIG. 4. Dependence of the fluctuations of $ on m/wipj. We set uj = 10 and e = 0.1. In (a) 
we show the evolution of r) for m/mpi = 0.5. The behavior of r) is quite different 
from that for small m. In (b) we draw A(t) for m/mpi = 0.1, 0.2, and 0.5. As m/rapi is 
larger, the fluctuations of $ become larger. 

FIG. 5. Dependence of the fluctuations of $ on the potential shape. We draw A(£) for 
e = 0.05, 0.1, and 0.5. We set uj = 10 and m = 0.1. Although the temporal behaviors of 
$(£, r) are quite different, the eventual behavior does not depend on the potential shape. 
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